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1.0  Introduction 


This  report  covers  the  period  from  January  1>  1963  to 
July  1,  1964.  A  bibliography  was  submitted  in  lieu  of  a  progress 
report  in  July  1963.  No  progress  report  was  submitted  in  January 
of  1964  since  several  major  pieces  of  work  were  very  near  com¬ 
pletion  and  technical  reports  have  been  prepared  on  these 
projects. 

J^The  theme  of  the  research  reported  here  1^  the  investigation 
of  basic  properties  of  semiconductor  materials  which  are  of  impor¬ 
tance  in  energy  conversion  applications,  and  the  study  of  physical 
mechanisms  involved  in  energy  conversion.  Much  of  our  recent 
effort  has  been  directed  toward  a  theoretical  understanding  of  a 
particular  s emic owttoc t or  (PbTe)  .  The  theoretical  calculations 
reported  here  have  suggested  several  experimental  investigations 
which  are  now  being  pursued^  The  theoretical  work  has  been  done 
in  conjunction  with  ProfesseV  Pratt's  Materials  Theory  group; 
this  collaboration  has  been  a  fruitful  one  and  we  expect  to 
continue  it  in  the  future.  « 

Section  2  of  this  report  deals  with  investigations  which 
have  been  completed.  These  are  described  fully  elsewhere  (in 
Technical  Reports  and  Theses)  and  the  abstracts  are  reproduced 
here  for  information. 

Section  3  presents  the  status  of  continuing  investigations, 
and  Section  4  gives  some  information  about  projects  which  have 
been  discontinued  as  unprofitable  at  this  time. 

A  list,  of  piiblications  for  the  period  covered  by  this  report 
is  attached. 
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2.0  Summary  of  Completed  Investigations 

The  research  reported  in  this  section  has  reached  a  stage 
where  a  formal  report  has  been  made.  In  most  cases,  this  Is  In 
the  form  of  a  distributed  technical  report;  occasionally,  the 
report  Is  In  the  form  of  a  thesis.  Abstracts  of  these  Inves¬ 
tigations  are  given  below. 

2.1  Mechanical  Energy  Flow  In  Crystal  Lattices,  L.  M.  Magld, 

Phys.  Rev.  134,  A158  (1964). 

(The  material  In  this  article  and  the  following  one  Is 
based  on  Technical  Report  No.  3,  May  12,  1962.  It  has 
been  rewritten  and  somewhat  extended.  Reprints  of  these 
two  articles  will  be  distributed  as  Technical  Reports 
when  available.) 

A  detailed  study  of  the  mechanical  energy  transported  through 
a  completely  general  crystal  lattice  In  the  classical  temperature 
range  Is  presented.  Both  a  mechanical  "Poyntlng"  vector  and  the 
corresponding  mechanical  *’Po)mtlng"  theorem  are  established  for 
an  arbitrary,  3 -dimensional,  anisotropic,  enharmonic  crystal 
lattice,  containing  many  atoms  per  cell,  Inclusive  of  substitu¬ 
tional  Impurities,  and  possessing  an  arbitrary  range  of  atomic 
Interactions.  The  extension  of  the  energy-flow  procedures  to  a 
quantum  treatment  Is  discussed,  and  the  resulting  quantum- 
mechanical  form  of  the  relation  of  mechanical  energy  flow  to  the 
group  velocity  for  the  restricted  case  of  a  perfectly  periodic 
and  harmonic,  but  otherwise  completely  general,  crystal  Is 
derived  for  the  classical  temperature  range. 

2.2  Attenuating  Wave  Analysis  of  Heat  Flow  In  Crystal  Lattices. 

L.  M.  Magld,  Phys.  Rev.  134,  A163  (1964). 

A  direct  approach  to  the  detailed  analysis  of  the  thermal 
relaxation  and  conduction  processes  In  crystal  lattices  In  the 
classical-temperature  range  Is  presented  In  terms  of  the  mechanical 
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energy  transported  by  attenuating  lattice  waves.  A  second-order 
classical  perturbation  procedure,  fonmilated  In  terms  of  tlme- 
and  space -dependent  normal  coordinates.  Is  used  to  solve  for  the 
dynamics  of  a  slightly  Imperfect,  nonlinear  general  crystal 
lattice  model  under  the  Influence  of  an  applied  temperature 
gradient.  Only  the  use  of  a  random-phase  assumption  for  Initial 
wave  amplitudes  at  t  ■>  0  and  statistical  averaging  of  the  subsequent 
dynamical  response  are  required  for  the  direct  determination  of 
.the  accepted  lattice  relaxation  times  from  the  time  dependence 
of  the  stored  mechanical -energy  density  (for  first-  and  second- 
order  perturbation  terms).  In  addition,  the  well-known  enharmonic, 
mass-fluctuation,  and  force-fluctuation  components  of  the  hlgh- 
temperature  thermal  conductivity  are  found  directly  from  the  steady- 
state  mechanical -power  density  within  the  lattice.  No  use  Is 
made  of  the  Boltzmann  transport  equation  or  standard  phonon 
scattering  theory,  although  the  results  obtained  are  wholly 
consistent  with  their  use.  Finally,  a  brief  discussion  Is  given 
on  the  extension  of  this  attenuatlng-wave  technique  to  the 
corresponding  quantum  treatment  of  low-temperature  heat  flow  In 
crystal  lattices. 

2.3  Relativistic  Effeots  In  Lead  Tellurlde,  J.  B.  Conklin,  Jr., 
Technical  Report  No.  8,  Nonr  1841(31),  June  1964.  (Also 
submitted  to  the  Department  of  Electrical  Engineering  In 
partial  fulfillment  of  the  requirements  for  the  degree  of 
Doctor  of  Science  at  the  Hassachusetts  Institute  of 
Technology.) 

The  Augmented  Plane  Wave  Method  for  calculating  the  band 
structure  of  a  solid,  has  been  extended  to  Include  the  relativistic 
terms  of  the  two -component  Hamiltonian  which  Is  obtained  by  two 
successive  applications  of  the  Foldy-Wouthuysen  transformation  to 
the  Dirac  Hamiltonian.  The  order  of  the  secular  equation  Is 
reduced  by  using  basis  functions  which  are  basis  partners  for  the 
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irreducible  representations  of  the  double  group.  These  functions 
are  constructed  from  the  solutions  of  the  non-relatlvlstlc 
Hamiltonian  for  lead  telluride.  The  forms  of  the  relativistic 
matrix  elements  between  them  are  found,  with  full  use  of  group- 
theoretical  techniques  to  simplify  the  calculation.  These  matrix 
elements  are  evaluated  numerically,  and  the  resulting  Hamiltonian 
matrix  Is  diagonalized  to  give  the  energy  levels  of  lead  telluride. 
The  resulting  band  structure  appears  to  be  reconcilable  with  the 
available  experimental  Information  about  lead  telluride  and.  It  Is 
hoped,  will  be  of  value  In  correctly  Interpreting  the  experimental 
data  for  the  material. 

2 .4  Calculation  of  Electronic  Properties  of  Strained  Lead 
Telluride .  L.  G.  Ferreira,  Technical  Report  No.  7, 

Nonr  1841(51),  April  1964.  (Also  submitted  to  the  Department 
of  Electrical  Engineering  In  partial  fulfillment  of  the 
requirements  for  the  degree  of  Doctor  of  Philosophy  at  the 
Massachusetts  Institute  of  Technology.) 

The  present  work  Is  a  study  of  the  electronic  properties  of 
Lead  Telluride  under  strain.  The  results  from  the  band  structure 
calculation  for  this  material  are  the  starting  point  for  a 
numerical  calculation  of  the  deformation  potentials  for  the 
Important  electron  states.  The  method  Is  based  on  Slater's  APW 
scheme.  An  expression  for  the  matrix  element  of  the  strain 
Hamiltonian  between  two  general  augmented  plane  waves  Is  found. 
Group  Theory  is  extensively  used  in  connection  with  the  reduction 
of  the  cubic  group  of  Lead  Telluride  by  a  uniaxial  strain,  and  In 
relating  different  deformation  potential  coefficients.  The  effects 
of  the  spln-orblt  Interaction  are  discussed  and  the  machine  program 
for  the  calculation  Is  explained. 

The  work  ends  by  comparing  the  calculated  deformation 
potentials  with  those  obtained  from  experimental  data  on  the 
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conduction  properties  of  the  material.  The  magnitude  of  the 
calculated  values  seems  to  confirm  the  order  of  the  energy  levels 
given  by  the  band  structure  calculation. 

2.5  Temperature  Dependence  of  the  Electron-Voltaic  Effect  In 

Silicon  P-N  Junctions,  Carl  E.  Rleck,  Jr.,  submitted  to  the 
Physics  Department  In  partial  fulfillment  of  the  requirements 
for  the  degree  of  Doctor  of  Philosophy  at  the  Massachusetts 
Institute  of  Technology. 

This  thesis  presents  a  theoretical  foundation  for  the 
electron -voltaic  effect,  beginning  with  the  mechanism  of  excita¬ 
tion  of  atomic  electrons  by  fast  beta  particles  and  ending  with 
the  temperature  dependence  of  the  behavior  of  the  carrier  pairs 
so  formed.  Experimental  results  are  also  presented  which  are 
substantially  In  agreement  with  the  theory.  From  the  experiments 
It  appears  that  correct  conclusions  are  reached  If  It  Is  assumed 
that  the  collections  of  electrons  and  holes  are  created  In  quasl- 
equlllbrlum,  even  thou^  It  seems  certain  on  other  grounds  that 
this  Is  not  the  case.  It  Is  concluded  that  quasi -equilibrium 
conditions  are  reached  In  times  which  are  short  In  comparison 
to  the  recombination  time. 

In  the  study  of  the  excitation  mechanism,  soft,  midrange, 
and  very  hard  collisions  are,  under  certain  conditions,  governed 
primarily  by  resonance  scattering.  The  relative  Importance  of 
this  effect  depends  mainly  on  the  off-diagonal  elements  of  the 
linear  momentum  operator.  Reliable  numerical  results  could  be 
obtained  only  If  these  elements  were  known.  For  this  reason  It 
is  suggested  that  the  electron-voltaic  effect  mi^t  be  a  useful 
tool  In  the  study  of  semiconductor  properties. 

In  the  treatment  of  carrier-pair  behavior  It  Is  predicted 
that  the  maximum  output  voltage  over  the  saturation  range  will 
be  nearly  linear,  decreasing  function  of  temperature.  On  the 
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basis  of  the  elementary  theory,  the  predicted  extrapolated  limit 
at  absolute  zero  is  the  gap  energy  divided  by  the  electron  charge, 
but  the  need  for  a  multiplicative  correction  Is  Indicated.  These 
predictions  are  confirmed  by  the  experimental  results. 

The  efficiency  of  an  electron-voltaic  cell  Is  limited  by 
the  Imperfect  utilization  of  bombarding  energy  In  carrier  pair 
production,  by  the  Imperfect  collection  of  carrier  pairs,  and  by 
a  factor  which  closely  resembles  a  Carnot  efficiency  at  low 
temperatures.  For  these  reasons  It  appears  that  the  electron- 
voltaic  effect  probably  cannot  be  made  to  provide  an  efficient 
means  of  converting  radioactive  source  power,  especially  at 
room  temperature. 

3.0  Current  Research 

3.1  Energy  Bands  and  Properties  of  Lead  Tellurlde 

The  calculated  electronic  energy  bands  reported  In 
Section  2.3  provide  a  model  which  may  be  useful  In  explaining 
various  experimentally  determined  properties  of  lead  tellurlde. 

There  Is  a  great  deal  of  experimental  Information  available  and 
It  should  be  possible  to  check  some  of  the  major  predictions  of 
the  model.  In  some  cases,  further  calculations  are  required  and 
In  others,  more  specific  experiments  need  to  be  carried  out. 

Some  of  the  possibilities  are  discussed  In  this  section. 

One  of  the  predictions  of  the  calculation  Is  the  size  of 
the  energy  gap  for  PbTe.  However,  as  pointed  out  In  Technical 
Report  No.  8,  the  gap  is  small  (^0.3  ev)  compared  to  the  total 
energies  calculated  8  ev)  and  the  calculation  is  not  sufficiently 
accurate  to  make  a  definite  prediction.  Instead,  the  energy  gap 
has  been  taken  as  a  parameter  to  be  adjusted  to  the  experimental 
value . 
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other  energy  differences >  which  could  be  detected  in  optical 
reflection  studies^  are  much  larger  and  one  would  expect  to  be 
able  to  check  the  theory  here.  However,  there  are  many  possible 
transitions  and  the  available  reflection  data  are  not  sufficiently 
precise  to  allow  an  unambiguous  interpretation. 

Using  the  calculated  energies  and  wavefunctions  (with  the 
energy  gap  adjusted  to  the  experimental  value),  Ferreira  and 
Pratt  have  calculated  the  effective  masses  for  the  valence  and 
conduction  bands  using  the  k*p  approximation.  They  obtain  very 
good  agreement  with  experimental  values.  Using  similar  techniques, 
they  have  calculated  g-factors  for  the  electrons  and  holes. 

These  factors  are  quite  different  from  the  free-electron  value 
and  it  may  be  possible  to  determine  the  values  experimentally. 

Such  experiments  are  being  considered. 

Another  check  on  the  theoretical  model  would  be  a  measurement 
of  the  work  function  or  electron  affinity  of  the  material.  With 
this  in  mind,  we  are  attempting  to  grow  epitaxial  films  of  PbTe 
and  perform  photoelectric  and  contact  potential  measurements  on 
these  films  without  removing  them  from  the  vacuum  chamber  in  which 
they  are  grown. 

One  may  obtain  some  information  about  the  symmetry  of  the 
wave  functions  associated  with  the  conduction  and  valence  band 
edges  by  looking  for  the  Knight  shift  in  nuclear  resonance. 
Published  data^  show  that  such  an  effect  exists  in  p-type  PbTe  and 
supports  the  theoretical  predictions.  Further  investigation 
using  this  technique  is  planned. 

In  order  to  use  the  model  to  interpret  transport  phenomena. 

It  will  be  necessary  to  study  the  scattering  processes  in  more 

^  I.  Weinberg  and  J.  Callaway,  Nuovo  Cimento  190  (1962). 
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detail.  In  particular,  we  are  investigating  the  electron -phonon 
interaction  theoretically.  Once  this  analysis  is  completed,  it 
should  be  possible  to  make  comparison  with  the  wealth  of  experi¬ 
mental  data  on  transport  processes. 

3.2  Heat  Flow  and  Attenuating  Wave  Analysis 

The  extension  of  the  analysis  reported  in  Sections  2.1  and 

2.2  to  the  low  temperature  region  where  quantum  effects  are 
expected  to  be  Important  is  nearing  completion.  The  results 
obtained  so  far  are  for  the  most  part  in  agreement  with  conven¬ 
tional  analyses,  but  the  present  method  gives  some  additional 
insight .  It  is  expected  that  a  complete  report  on  this  work 
will  be  published  in  the  near  future. 

4.0  Discontinued  Investigations 
4.1  Energy  Bands  in  Bl2Te^ 

We  have  considered  a  variety  of  experimental  techniques 
which  might  be  useful  in  determining  in  more  detail  the  energy 
bands  of  Bi2Te2.  In  particular,  we  were  interested  in  the 
departure  of  the  bands  from  parabolic  in  state  density,  and  in 
the  reasons  for  the  temperature  dependence  of  the  Hall  effect . 

After  a  detailed  consideration  of  magneto -optical  effects, 
it  seems  that  the  Interpretation  of  these  experiments  is  not 
sufficiently  precise  for  us  to  obtain  the  desired  information. 
Part  of  the  problem  lies  in  the  nature  of  the  material  and  the 
available  samples;  using  realistic  estimates  for  these  samples, 
it  was  clear  that  the  desired  information  would  be  obscured  by 
Other  effects.  In  some  cases,  the  experiments  were  so  difficult 
that  it  does  not  seem  worthwhile  to  attempt  them.  We  have  there¬ 
fore  suspended  our  efforts  in  this  direction. 
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4 . 2  Effects  of  Strain  in  PbTe 


We  have  carried  out  a  preliminary  investigation  of  the  effects 
of  strain  on  PbTe  tunnel  diodes.  Several  PbTe  tunnel  diodes  were 
constructed  and  measurements  of  the  effects  of  strain  were  made 
on  the  best  of  these.  The  diodes  were  made  on  material  oriented 
so  that  the  largest  effects  of  strain  were  anticipated. 

These  diodes  do  not  have  characteristics  which  follow  any 
of  the  present  theories  of  tunneling  current.  At  best,  if  the 
effects  of  strain  were  large,  it  might  be  possible  to  relate  the 
strain  qualitatively  to  a  parameter  of  the  diode.  However,  the 
effect  of  strain  was  a  small  one  and  did  not  occur  in  a  region  of 
the  characteristic  which  could  be  easily  interpreted.  Since  our 
major  interest  is  in  the  material  and  not  in  the  theories  of 
tunneling,  we  have  abandoned  this  particular  line  of  investi¬ 
gation  of  PbTe. 
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